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Q. 

^ . Quantum fluctuations of electromagnetic radiation pressure are discussed. 

£S) | We use an approach based on the quantum stress tensor to calculate the 

fluctuations in velocity and position of a mirror subjected to electromagnetic 
radiation. Our approach reveals that radiation pressure fluctuations are due 
to a cross term between vacuum and state dependent terms in a stress ten- 
sor operator product. Thus observation of these fluctuations would entail 



CN ! experimental confirmation of this cross term. We first analyze the pressure 



o 



fluctuations on a single, perfectly reflecting mirror, and then study the case 
of an interferometer. This involves a study of the effects of multiple bounces 
in one arm, as well as the correlations of the pressure fluctuations between 
arms of the interferometer. In all cases, our results are consistent with those 
previously obtained by Caves using different mehods. 

I. INTRODUCTION 

> : 

^ . Classically, a beam of light falling on a mirror exerts a force and the force can be written as 

the integral of the Maxwell stress tensor. When we treat this problem quantum mechanically, 
then the force undergoes fluctuations. This is a necessary consequence of the fact that 
physically realizable quantum states are not eigenstates of the stress tensor operator. These 
radiation pressure fluctuations play an important role in limiting the sensitivity of laser 
interferometer detectors of gravitational radiation, as was first analyzed by Caves His 
approach is based on the photon number fluctuations in a coherent state, and we will refer 
to it as the photon number approach. The purpose of this paper is to examine radiation 
pressure fluctuations using the quantum stress tensor. This requires the correlation function 
of a pair of stress tensor operators, which will be discussed in Sect. P|. There it is shown 
that the correlation function can be decomposed into three parts, a term which is fully 
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normal-ordered, a state independent vacuum term, and a "cross term" which can be viewed 
as an interference term between the vacuum fluctuations and the matter content of the 
quantum state. It is this cross term which will be of greatest interest in this paper, as it 
is responsible for the radiation pressure fluctuations in a coherent state. These fluctuations 
will be discussed in Sect. [TTI| for the case of a laser beam impinging upon a single, perfectly 
reflecting mirror. The analysis will be done first using the photon number approach, and 
then using the stress tensor approach. In the latter case, we show how the calculations may 
be performed in coordinate space, where an integration over space and time is needed to 
remove a singularity in the cross term. We then show how to obtain the same result more 
simply using a othonormal basis of wavepacket modes. In Sect. [Ill D , we examine the case of 
a single mode number eigenstate, and show that the radiation pressure fluctuations vanish. 
In the stress tensor approach, this arises from a cancellation between the cross term and 
the fully normal ordered term. In Sect. [V], we turn to the discussion of an interferometer. 
We first study the effects of mutiple bounces in a single arm, and reproduce the result 
that the effect of the radiation pressure fluctuations grows as the square of the number of 
bounces. We also show how our approach may be used to discuss the situation where the 
interferometer arms are Fabry-Perot cavities. Finally, we discuss the correlation between 
the fluctuations in the two interferometer arms, and show why they are in fact uncorrelated. 
Our results are summarized and discussed in Sect. |V[ 



II. ENERGY-MOMENTUM TENSOR FLUCTUATIONS 

It is well known that stress tensor operators can be renormalized by normal ordering: 

: Tiiv '■= T^ v — (T M ;,)o , (1) 

which is subtraction of the Minkowski vacuum expectation value. However, the quantity 
(: T^ u (x) :: T pa (x') :) is still divergent in the limit that x' — > x. The divergent part of this 
quantity can be decomposed into a state-independent part and a state-dependent part. To 
do so, we may use the following identity, which follows from Wick's theorem 

: X 2 :: 3 04 : = : 01020304 : 

+ : 0103 : (0 2 04)o+ : 0104 : (0203)o 
+ : 0203 : (0i04)o+ : 0204 : {4>i4>z)o 

+ (0103)o(0204)o + (0104)o(0203)o • (2) 

Here the 0j are free bosonic fields and ( )o denotes the Minkowski vacuum expectation value. 
The first term is fully normal-ordered, the next four are cross terms and the final two are 
pure vacuum terms. The physics of these various terms was discussed in Ref . || . Here 0i and 
02 are evaluated at point x, whereas 03 and 04 are evaluated at point x'. In the coincidence 
limit, x' — > x, the fully normal-ordered term is finite, but the cross term and vacuum terms 
diverge. The singularity of the cross term is of particular significance because, unlike the 
vacuum term, it is state-dependent. The fully normal-ordered term will not contribute to 
the fluctuations so long as the quantum state is a coherent state. The pure vacuum term 
will also not contribute so long as we restrict our attention to the differences between a given 
quantum state and the vacuum state. 
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If the quantum state is other than a coherent state, there are also state-dependent stress 
tensor fluctuations in the fully normal ordered term. These fluctuations were discussed in 
Refs. [|]-|7j], especially in a context where the stress tensor is the source of gravity. The 
normal ordered term is always finite and does not present a divergence problem, in contrast 
to the cross term. The latter term can only be made meaningful if one examines space or 
time integrated quantities and has a prescription for defining the resulting integrals. We 
may schematically express the expectation value of a product of stress tensor operators as 

(• Tpi/ .. Tp a .) (. TpJTpa •) ~\~ ^TpyTpcr) cross {Tpi/Tp a ^ Q , (3) 

where the three terms of the right hand side are, respectively, the fully normal ordered term, 
the cross term, and the vacuum term. For a single mode coherent state \z), 

(z\ : Tp V T pa : \z) = (z\ : Tp U : \z)(z\ : T pu : \z) . (4) 

In such a state, the fluctuations of the stress tensor are described by quantities of the form 

(AT 2 ) = (: Tp U : 2 ) - (: T, v :) 2 = (: T, v : 2 ) cross + (: T, u : 2 ) . (5) 

If we are interested only in the changes in (AT 2 ) when the quantum state is varied, then 
the pure vacuum term can be ignored, and only the cross term is important: 

(AT 2 )^(:T^: 2 ) cross . (6) 

Note that we do not mean to suggest that there is no physical meaning to the pure 
vacuum term. It presumably describes fluctuations of the stress tensor components in the 
Minkowski vacuum state. More precisely, if one measures a spacetime averaged component, 
the result of the measurement should undergo fluctuations which vary as an inverse power 
of the size of the averaging region. However, in a non-vacuum state, the magnitude of the 
cross term will grow as the mean energy density in the state. Thus there will be a regime 
in which the effects of the cross term dominate those of the vacuum term. 

We must resolve the issue of the state-dependent divergences in the cross term if it is to 
have any physical content. This issue was discussed by us in Ref. ||, where it was shown that 
although the stress tensor correlation function is singular in the coincidence limit, integrals 
of this function over space and time can still be well-defined. The cross term in the stress 
tensor correlation function has the form 

(Tpu{x)Tp a (x ))cross — J~ ~ a/)^ ' ^ 

where F(x,x') is a regular function of the spacetime points x and x', and (x — x') 2 is the 
squared geodesic distance between them. Integrals of the correlation function appear to be 
formally divergent, but nonetheless may be defined by an integration by parts procedure. 
Suppose for the sake of illustration that the integrations are over time only, and note that 

ln(t - t'f . (8) 



(t - t'Y 12 dt 2 dt' 2 

If the function F vanishes sufficiently rapidly at the endpoints of the integrations, we can 
write 
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( dtdt' ^^'l =-— [ dtdt'ln(t-f) 2 ^--^-F(t,t') . (9) 
7 (t - t') 4 12 7 K J dt 2 dt' 2 K J K ' 

This procedure provides a way to define integrals with singular integrands, and has been 
discussed by various authors [@|§. As we will see below in Sect. [111 B| and in the Appendix, 



the integrals of the cross term which describe radiation pressure fluctuations can be made 
finite by a similar procedure. 



III. INDUCED MOMENTUM FLUCTUATIONS OF A SINGLE MIRROR 

It is well known from classical physics that a beam of light falling on a reflecting or 
absorbing surface exerts a pressure. This pressure may be computed by integration of 
the appropriate component of the Maxwell stress tensor over the surface. It may also be 
computed by counting photon momenta. Let us illustrate the latter method, which we 
will call the "photon number" approach. If an incident monochromatic beam of angular 
frequency uj and energy density p strikes a surface, the mean number of photons striking 
per unit time per unit area is[] p/u. If the light is perfectly reflected, each photon imparts 
a momentum 2uj to the surface, resulting in a radiation pressure of 2p. As expected, both 
the stress tensor and the photon number approaches yield the same answer. 

However, these calculation give only a mean value. The radiation pressure should undergo 
fluctuations about this mean. In the photon number viewpoint, these fluctuations arise from 
fluctuations in the rate of photons striking the surface. In the stress tensor viewpoint, the 
fluctuations arise because the quantum state of the radiation field is not an eigenstate of 
pressure. The main purpose of this section is to examine radiation pressure fluctuations 
in a single mode coherent state from both viewpoints, and to compare the results. In 
one subsection, Sect. |IIID| , we will also examine the case of a single mode photon number 
eigenstate. 



A. The Photon Number Approach 

In this approach, the radiation pressure fluctuates because of statistical fluctuations 
in the numbers of photons striking the surface. Suppose that a beam of light with angular 
frequency uj is described by a single mode coherent state, \z), an eigenstate of the annihilation 
operator, a\z) = z\z). The mean number of photons which strike a mirror in time r is 

( n ) = ( a t a ) = \z\ 2 . (10) 

If the mirror is perfectly reflecting, then the mean momentum transferred is the expectation 
value of the operator 

p = 2un . (11) 



Units in which K = c = 1 will be used throughout this paper. Electromagnetic quantities are in 
Lorentz-Heaviside units. 
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The dispersion of this momentum is given by 

(Ap 2 ) = (p 2 ) - (p) 2 = Aco 2 ((n 2 ) - (n) 2 ) . (12) 

In a coherent state, 

(n 2 ) - (n) 2 = (afoota) - {a ] a) 2 = ((a^V) + (a f a) - (a f a) 2 = (a^a) = (n) . (13) 

Thus 

(Ap 2 ) =4uj 2 (n) =4luA P t, (14) 

where p is the mean energy density of the incident beam, and A is its cross sectional area. 
If the mirror is a free body with mass m, the mean squared velocity fluctuation is 

(Av 2 )=4^fr. (15) 



B. The Stress Tensor Approach 

An alternative approach to the problem of radiation pressure fluctuations is the method 
of stress tensor fluctuations. It is well known that one can calculate the force on a surface 
by integration of the relevant component of the stress tensor over that surface. It thus 
seems reasonable to expect that the fluctuations in this force can also be computed from 
the quantum stress tensor. There is, however, a problem which needs to be resolved in this 
approach. This is that products of stress tensor operators are not well defined at coincident 
points. Even the integrals of these products are formally divergent and need a regularization. 
A regularization scheme was used in our previous work || and will be adopted again in this 
paper. The main idea is to treat the measuring process as a switch-on-switch-off process 
and do an integrations by parts. 

Consider a mirror of mass m which is oriented perpendicularly to the x-direction. If the 
mirror is at rest at time t = 0, then at time t = r its velocity in the x-direction is given 
classically by 

v = — / dt [ daT xx , (16) 



m Jo J a 



where is the Maxwell stress tensor, and J A da denotes an integration over the surface of 
the mirror. Here we assume that there is radiation present on one side of the mirror only. 



Otherwise, Eq. (16) would involve a difference in T xx across the mirror. When the radiation 
field is quantized, is replaced by the normal ordered operator : :, and Eq. ( JTB| ) becomes 
a Langevin equation. The dispersion in the mirror's velocity becomes 

(Av 2 ) = i- f dt { T dt' ( da ( da'[(: T xx (x) :: T xx {x') :) - (: T xx (x) :)(: T xx (x') :)] (17) 



m? Jo Jo J a J A 



As discussed above, when the quantum state of the radiation field is a coherent state and 
we ignore the pure vacuum term, then the dispersion in T xx is given by the cross term alone, 
and 
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(Av 2 ) = — \ T dt f dt' ( da ( da'(T xx (x)T xx (x')) cross . (18) 
m z Jo Jo J a J A 

The components of the energy-momentum tensor for the electromagnetic field are (Lorentz- 
Heaviside units are used here.) 

T ij = S ij T oo _ ( E i E j + B i B jj ^ ^g) 

T° l = e ijk E j B k , (20) 

and 

T 00 = ^{E 2 + B 2 ) . (21) 

Here E l and B l are Cartesian components of the electric and magnetic fields, respectively. 
In particular, 

T xx = l(E 2 + E 2 + B 2 + B 2 z ). (22) 
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We now assume that a linearly polarized plane wave is normally incident and is perfectly 
reflected by the mirror. Take the polarization vector to be in the y-direction, so that 
E z — By — 0. At the location of the mirror, E y — 0, and only B z contributes to the 
stress tensor. Thus, when we apply Eq. (0) to find {T xx (x)T xx (x')) cross , the only nonzero 
quadratic normal-ordered product will be (: B z (x)B z (x') :). The result is 

(T xx (x)T xx (x')) cross = (: B z (x)B z (x') :)(B z (x)B z (x')) . (23) 

The ^-component of the magnetic field operator may be expressed in terms of mode functions 

as 

B z (x) = '£(a u B u + atB* u ), (24) 

where the mode function is 

B (J = -2Cca&(ux)e-** t . (25) 
Here C is the coefficient for box normalization in a volume V 



C = J£. (26) 
The coherent state is an eigenstate of the annihilation operator 

a u r\z) = 5 w/LlJ z\z) , (27) 

where z is a complex number 

z=\z\e-*+. (28) 
The expectation value of the normal ordered product of field operators is now 
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(z\ : B z (ti)B z (t 2 ) : \z) = 1QC \z\ cos(u;ti + <p) cos(ut 2 + 4>) cos(a>a;i) cos(a;a;2) (29) 

The vacuum magnetic field two-point function in the presence of a perfectly reflecting plane 
at z = is given by 

{B z {t±, Xi)B s (t 2 , x 2 )) = {B z {t u yn)B z {t 2 , x 2 )) S0 + {B z {t u x x )B z {t 2 , x 2 )) /0 . (30) 

The first term is the two-point function for empty space, 

m(f v\R(i - (^i-^) 2 + |xi-x 2 | 2 -2(2: 1 -^) 2 

{B z (tx, Xi)B z (t 2 , x 2 J) B0 - — — j — . (61) 

The second term is an image term 



(S,(ti,x 1 )B z (t 2 ,x 2 ))jQ = (B x (t 1 ,x. 1 )B s (t 2 ,x. 2 )) E0 



(32) 



Both terms give equal contributions to the radiation pressure fluctuations on a mirror located 
at z — 0. 

We can now combine these results to write Eq. (ITS) as 



2\ 32C^ 1^1 f i f i / \ 

^ v ) — — o — I dai / da 2 J (33) 

TT in J A J 



where 



with 



and 



J — dti dt 2 w . 1 - l„ — r^r cos(utt) cos(ujt 2 ) , (34) 
Jo Jo l(ti — t 2 y — b l \ A 



zi - z 2 f - ( Vl - y 2 f (35) 



b 2 = (yi-y 2 ) 2 + (zi-z 2 ) 2 . (36) 

Here we have set = arg(^) = 0, as it just shifts the origin of time, and taken the location 
of the mirror to be at x = 0. The integral J is evaluated in the Appendix in the limit of 
large r, with the result 

J = { \b 2 (b 2 + a)u 2 + b 2 - 3a] sinbu + ub{3a - b 2 ) cosku} . (37) 

The singularities in the integrand of Eq. (|34T ) are third order poles, which are evaluated 
using an integration by parts. 

We next need to perform the spatial integration over the area of the mirror which is 
illuminated by the laser beam. Assume that the illuminated region is a disk of radius R and 
hence area A = irR 2 , and that the incident flux is uniform over this disk. (This assumption 
is not essential, but simplifies the calculations.) If we take the origin for the da\ integration 
to be at y 2 = z 2 — 0, then 
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1= da 2 dai [b 2 (b 2 + a)cu 2 + b 2 — 3a] smbcu + ujb(3a — b 2 ) cosbu] b 

J A J A 



R 



2- 



r dr I d6 



dd2 

ur{r 2 — 3r 2 sin 2 6) cos(a;r) + [3r 2 sin 2 6 + r A u 2 — r 2 (l + r 2 sin 2 6uj 2 )] sin(u;r) 
R (1 + u 2 r 2 ) sinfcur) — ur cosiur) 



7T da 



dr . 



(38) 



Further assume that R ^> uj 1 . Then we can let R — > oo in the upper limit of the r 
integration. The da% integration now simply contributes a faxtor of A, and we have 

(1 + uj 2 r 2 ) sin(a;r) — ur cos(a;r) 



TlUlA 



dr 



[1 + it 2 ) sin(-u) 



MCOS (U) 



du 



du 



d ( sin u 



smu — 



du\u 



= txujA 
= ItxujA . 

In the last step, we used J °° du sinw = lim^o /o°° du sinwe - "" = 1 . Thus we have 

1QC 2 \z\ 2 uAt 



The energy density in the incident wave can be written as 

12 



P 



10\Z\ 



2C 2 \z 



2UI2 



so we can express the velocity fluctuations as 



Note that this result agrees with that from the photon counting approach, Eq. 



(39) 



(40) 



(41) 



(42) 



C. The Wavepacket Approach 

Here we wish to provide an alternative derivation of the momentum fluctuations of a sin- 
gle mirror using the stress tensor approach. Rather than performing all of the calculations in 
coordinate space, as was done in Sect. [Ill rj| , we will use an approach based upon wavepacket 
modes. This approach will prove useful in discussing interferometer noise. Assume that the 
occupied mode is a wavepacket which is sharply peaked at frequency u. Using Eq. (|TJ) and 
Eq. ( p3|) for a coherent state, the momentum fluctuation becomes 

(Ap 2 ) =m 2 (Av 2 } 

= J dtda dt> da' (T xx (x)T xx (x')) cross 

= J dtda J dt'da'(: B z (x)B z (x') :) (B z (x)B z (x')) . (43) 
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Let the magnetic field operator B z (x) be expanded in terms of a complete set of positive 
frequency wavepacket modes {uj(x)} : 

B z {x) = ^{aj Uj(x) + OjU*(x)] . 
j 



(44) 



For our purpose, we take these modes to be fairly sharply peaked in frequency and use the 
normalization condition 



u)u y d 3 X = -Uj Sjj, , 



(45) 



where utj is the mean frequency of packet j. More generally, we should expand the vector 
potential A in terms of wavepackets jf 

A = YJv&h + H - c -) > 

j 

where (fj, fji) = djji and (fj, fj') is the Klein-Gordon inner product. Then the modes Uj are 
expressed as derivatives of the modes fj. Consider a single mode coherent state \z) as the 
quantum state, and let the mode be a wavepacket Uq. Then with a suitable choice of the 
phase of this mode function, we can write 

(: B z (x)B z (x') :) = \z\ 2 [u (x) + u*(x)}[u (x') + u* (x')] . (46) 

Note that the integrals in (Ap 2 ) are of the form / dtda rather than / d 3 x . If the integrand 
is a function of t — x alone, or t + x alone, these are equivalent: 



/oo /*oo /*oo 

duf{u)= I duf(—u)= / dxf(t — x) 
-co J —oo J —CO 



and 



/OO /"OO 
dtf(t + x) = / dxf(t + 
-OO J —CO 



X) 



(47) 



(48) 



where u = t — x. However, when the mirror is present, Uq contains pieces moving in both 
directions: 



u = u 0I + u 0Ri 



(49) 



where uqi is the incident wavepacket and uqr is the reflected wavepacket. The key feature 
that we will use is that uq is orthogonal to Uj (j ^ 0), but uqi and uqr are not orthogonal 
to each other. Inserting Eq. (|44]) and Eq. (f46D into Eq. ([43D yields 

(Ap 2 ) = \z\ 2 dtda / dt' da'[uo(x) + ul(x)][uo(x f ) + u* (x')] 



- Y\uj{x)u*{x') + U*(x)Uj(x')\ 

z j 

:\ 2 ^2Re J dtda J dt' da f u (x)uQ(x')u*(x)uj(x') 
j 

:\ 2 ^2Re J uo(x)u*(x)dt da' J ul(x')uj(x')dt' da 



uo(x)u* Q (x)dt da 



1 2 



(50) 
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Using Eq. fl49|) , the integral in the last expression becomes 

J = y u (x)ul(x)dt da 
= J dtda (u 0I + u b)(u* oi + Uq R ) 

= J dtda (\u i\ 2 + \u 0R \ 2 + u 0I ul R + uIjUqr) . (51) 

Here the only difference between u i and u Q r is their direction of travel. If u j = f(t — x), 
then Uqr = f(t + x), and Uqj = Uqr at the mirror, and 

J dtdauQ[U* QR = J dtdau* 0I UQR = J dt da \u i\ 2 = J dt da \u r\ 2 . (52) 

Thus 

I — J u (x)ul(x)dt da = A J dtda\uoi\ 2 = A J d 3 x\uoi\ 2 = 2lu , 

and the momentum fluctuation becomes 

(Ap 2 ) =Aoo 2 \z\ 2 = Aoo 2 (n). (53) 

This is same as the result of the photon number counting approach, Eqs. (pT5[), and of the 
integration by parts in Sect. |l II B|, (|^). 



D. A Number Eigenstate 

Most of this paper deals with single mode coherent states. However, in this subsection, 
we wish to turn aside from the main line of development and discuss the case of a single 
mode number eigenstate. It is apparent in the photon number approach that there should 
not be any radiation pressure fluctuations in such a state. In the stress tensor approach, the 
situation is less clear, as both the fully normal-ordered term and the cross term are nonzero. 

We assume that the quantum state \n) is a number eigenstate of a single mode. Expand 
the magnetic field operator in terms of complete set of mode functions using Eq. (^) to find 
the expectation value of the fully normal ordered term 

(n\ : T xx (x)T xx {x') : \n) 
= \(n\:B 2 z {x)B 2 z {x'):\n) 

= \<n - l)[B* 2 {x)Bl{x') + Bl{x)B* 2 {x>) + 4|^(o;)| 2 |i?.(a;')| 2 ] . (54) 



Here B u is the mode function for the occupied mode, and is assumed to be given by Eqs. ( |25|) 
and (|26|) . A similar procedure leads to the result for the cross term 

. T xx (x} .. T xx {x ) . cross 

= (n\ : B z (x)B z {x') : \n) (B z (x)B z (x')) 

= n [B u {x)Bl{x') + B*(x)B u {a/)] (B z (x)B z (x')) . (55) 
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n(n — I) 



The mean pressure is 

(: T xx (x) :)„ = (n\ : T xx (x) : \n) = n\B u (x)\ 2 . (56) 

The momentum deviation due to the fully normal ordered term becomes 

(n| : Ap 2 : \n) n = (: p 2 :) n - (: p :) 2 n 

dtda I dt'da' ({: T xx (x)T xx (x') :) n - (: T xx (x) :) n (: T xx (x') :)„) 

dtda J dt'da' (B*J(x)Bl(x') + Bl(x)B*J(x')) 
-n J dtda J dt! da' \B w [x)f\B u {x')f . (57) 

Integrals such as / dt B 2 (x) contain rapidly oscillating integrands. We will assume that these 
integrals average to zero and can be ignored. The remaining integrals, involving l-B^I 2 are 
straight forward. The result for the contribution of the fully normal ordered term to the 
momentum deviation is 

(: Ap 2 :)„ = -Anu 2 . 

The cross term is 

(n\Ap 2 \n) cross = J dtda J dt' da' (n\ : B z {x)B z {x) : \n)(B z (x)B z (x'))o 

= Ancu 2 . (58) 

Note that the cross term calculation is almost identical to the calculation in the case of a 
coherent state. The momentum deviation due to the fully normal ordered term and the 
cross term cancel out each other and yield zero momentum deviation 

(n| Ap 2 \n) = (n\ : Ap 2 : |n) + (n| Ap 2 \n) cross = , (59) 

which is expected in photon number approach. This calculation shows the agreement be- 
tween the photon number and the stress tensor approaches. In the stress tensor approach, 
the fluctuations in the normal ordered term are anticorrelated with those described by the 
cross term. 



IV. NOISE IN AN INTERFEROMETER 



A primary application of the result in Sect. [Ill A] is to estimate the radiation pressure 
noise in an interferometer. (See Fig. |l|) The laser beam bounces b times in a arm of the 
interferometer before being recombined. The masses at the end of each arm are subject to 
velocity and position uncertainty due to the radiation pressure fluctuations. 
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FIGURES 




FIG. 1. The essentials of a Michelson interferometer. An input laser beam ("In") is split by the 
beam splitter S. The split beam bounces b times (here b = 3) in each arm between a free mass (mi 
or 7712) and a fixed mirror (Mi or M2) before being recombined to form the output beam ("Out"). 
The mirrors of interest for the radiation pressure fluctuations are located on the free masses. 

A. Position Uncertainty in the Photon Number Approach 

In this subsection, we will review the conventional, photon number approach to calculat- 
ing interferometer noise . The effect of multiple bounces is accounted for by multiplying 
the momentum operator in Eq. ( p]) by a factor Of b 



p = 2bun . 

This introduces a factor of b 2 in the mean squared velocity fluctuation 



<A, 2 )=4^r. 

m 2 



(60) 



(61) 



The root mean squared position uncertainty of each mirror due to radiation pressure fluc- 
tuations is then of order 



Ax 



rp 



b T2 



m 



(62) 



where P = Ap is the mean power in the laser beam. There is another source of noise in 
the interferometer, the photon counting error, also known as shot noise. This arises from 
the uncertainty in the location of an interference fringe when a finite number of photons are 
counted, and is of order 



Ax 



pc 



2bVujPr 

If we minimize the net squared position uncertainty 



(63) 



12 



{Ax rp f + {Ax pc f (64) 
with respect to P, the result is the optimum power 



m 



opt 2ur 2 b 2 

When P = P op t, the position uncertainty becomes 



(65) 



Ax = Ax sq i = J— , (66) 
V m 

known as the standard quantum limit. 

The standard quantum limit is the position uncertainty which one obtains after time r 
by preparing a particle in a wavepacket state with initial spatial width Axo and momentum 
spread Ap = mAxo/r. After time r has elapsed, the spread in the width of the wavepacket, 
Apr/m, is of the order of Axo- Thus the standard quantum limit can be interpreted as 
being the minimum position uncertainty which can be maintained for a time of the order of 

T. 



B. Position Uncertainty in the Stress Tensor Approach 

Here we wish to look in more detail at how the position uncertainty arises in a descrip- 
tion involving the quantum stress tensor. The position of each mass is disturbed by the 
pressure fluctuations, leading to an error in the measurement of the optical path length. 
The displacement of a given mass is described by the time integral of Eq. ( |16|) 



x — x(0) = — / / / T xx (xi)daidt\dt (67) 
m Jo Jo J a 



The dispersion in position is then given by an expression analogous to Eq. (|T 

(Ax 2 ) = (x 2 ) - (x) 2 

= -K f [ da x da 2 [ [ dtdt' [ [ [(: T xx {x x ) :: T xx (x 2 ) :) 
m 1 J A J A Jo Jo Jo Jo 

- (: T xx ( Xl ) :) (: T xx (x 2 ) :)] dt x dt 2 . (68) 

As before, only the cross term in the stress tensor fluctuations will contribute when the 
quantum state is a coherent state. If we use this fact, and then take the second derivative 
with respect to r, we can write 



d 2 2 
— (Ax 2 ) = — 
dr 2 ' m 2 



dh J dt 2 J d^ J da 2 (T xx (ti)T xx (t 2 )} ( 

- / dt' [ dt 2 [ dai [ da 2 (T xx (T)T xx (t 2 )) crost 
o Jo J J 



(69) 



Now we need to assume that the laser beam is switched on in the past and then switched 
off in the future. This issue was discussed in Ref. ||, where integrations by parts were 
performed in order to deal with the singular behavior of the cross term. The asymptotic 



13 



condition insures that the surface terms arising in the integrations all vanishes. In the present 
calculation, we require that the normal ordered factors vanish at t = r, and hence the second 
term on the right hand side of Eq. fl69|) vanishes. The remaining term is proportional to 
(A-u 2 ), so that 

^(Ax 2 ) = 2(A^ 2 ). (70) 

Thus if (Av 2 ) oc r, then (Ax 2 ) oc r 3 . This calculation is the justiication for Eq. (p2]). 

There are two major issues to be studied in remainder of this section. One is to study 
the effects of multiple bounces within one arm of the interferometer, which will be done in 
the following subsection. The other is to study whether there are any correlations between 
the two arms. The discussion in Sect. [IV A| implicitly assumed the absence of correlations, 



and a justication for this assumption was given by Caves 0. In the present context, there is 
a correlation term (Ax\ 2 ) which takes the form of Eq. (|68|) with X\ on one mirror and x 2 on 
the other mirror. We will argue below in Sect. IV D 2| that this term is negligible compared 
to (Ax 2 ) for each arm separately. 



C. Multiple Bounces in One Arm 

1. A Delay Line 

Now we wish to consider the situation where a laser beam bounces several times between 
a pair of mirrors, as illustrated in Fig. 1. This arrangement is sometimes called a "delay 
line". Suppose that the beam is recycled b times within a single interferometer arm. We 
have already seen that in the photon number approach, the momentum fluctuation of the 
end mirror, (Ap 2 ) is now proportional to b 2 . Specifically, 

(Ap 2 ) = b 2 (Ap 2 ) u (71) 

where (Ap 2 )i is the single bounce result given in Eq. (|53|) . One can understand this result in 
the following way: if there is a fluctuation in the number of photon entering the interferom- 
eter arm, that fluctuation is maintained on each of the successive bounces. If slightly more 
than the expected number of photons hit the mirror on the first bounce, the same excess 
will reappear on later bounces. One can picture the same photons as simply recycling b 
times. However, in the stress tensor approach, it is less obvious how the b 2 factor will arise. 
This factor requires that the stress tensor fluctuations at the different spots on the mirror 
be exactly correlated with one another. 

We can understand this by returning to Eq. (|50"D. If the integrations run over all time and 
over the entire area of the mirror, then they pick up contributions from all of the bounces. 
Let us first suppose that the spots formed on succesive bounces overlap on the same region 
of the mirror. In this case, the mode function uq is approximately periodic for b periods: 

Uo(t,x) ^u (t + 2Ln,x.), n = 1, 2, • • • , b - 1 . (72) 

Let the first bounce occur at time t — 0, and subsequent bounces at t — 2Ln. More 
precisely, these are the mean times at which the wavepacket hits the mirror. Let T be some 
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time interval which is long compared to the length of the wavepacket, but short compared 
to 2L. Equation (51) becomes 



(Ap 2 ) = \z\ 2 l\ (73) 

where 

1=1 dt I da u (x)uq(x) = V] / dt dauQ(x)u* {x) . (74) 

J-oo J n=0 J2Ln~T J 

Note that the time intervals which are ignored in going from the first form to the second are 
ones in which uq = on the mirror, that is, in between bounces. However, the periodicity 
property, Eq. (J72J), implies that each term in the sum is equal. Furthermore, each term gives 
the same contribution to / as was found in the single bounce case: 

I = bJ dt J dau (x)ul(x) = 2bu . (75) 

Thus we obtain Eq. (|7"T|). 

Note that it does not matter whether the spots formed on the various bounces actually 
overlap on the mirror or not. If they do not, then Eq. ( |72|) is replaced by a more complicated 
relation involving an offset in position for the different bounces. However, once the area 
integration is performed, this is irrelevant, and we still obtain Eq. (0). Note that we are 
assuming that on all bounces, the beam is nearly perpendiular to the mirror. 



2. A Fabry-Perot Cavity 



u, 



Ru, 



R 2 u 



M m 

FIG. 2. A Fabry-Perot cavity. Here the right mirror m is perfectly reflecting, and the left 
mirror M is partially reflecting, with a reflection amplitude R. An initial wavepacket uo returns 
to m as Ruq on the first bounce, as R 2 uq on the second bounce, ect. 
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The delay line arrangement sketched in Fig. [I] and discussed above implies a precise num- 
ber of bounces. Another possibility, which is more likely to be used in actual interferometers, 
is the Fabry-Perot cavity, illustrated in Fig. |2|. Here at least one of the mirrors is partially 
reflecting, leading to a finite storage time for a wavepacket in the cavity. We will discuss the 
case where the mirror on the free mass is assumed to be perfect, but the opposite mirror 
in the cavity is not. This assumption allows us to continue to use our previous expressions, 
especially Eq. (ffij). If the mirror on the free mass is not perfect, then it is necessary to 
modify this expression and include electric field terms as well. 

Let R be the complex reflection amplitude for the imperfect mirror, so that \R\ 2 is the 
fraction of the power reflected on each bounce. We assume that once inside the cavity, a given 
wavepacket mode bounces an infinite number of times, but with diminishing amplitude. The 
effective number of bounces, b', can be defined by 

\R\ 2b ' ' = -• (76) 

Thus the energy stored in an occupied wavepacket mode is reduced by a factor of \ after b' 
bounces. We can then write 

|i?| 2 = e-^^l-^, (77) 

The effect of the finite reflectivity of the left mirror is to introduce a factor of r each 
time the wavepacket returns to the right mirror. Recall that the magnetic field mode Mo has 
no phase shift upon reflection from the perfect mirror. We can express this as the following 
condition on the mode function: 

u (t + 2Ln,x) = R n u (t,x), n = 0, 1,2,---, -T < t < T . (78) 

That is, «o(t,x) for —T < t < T is the initial form of the wavepacket when it hits the left 
mirror for the first time. The above relation gives its form when it returns for the n-th time. 
We can now use this relation to write the analog of Eq. ( |7^ ) : 

I = dt dau (x)u* (x) = V / dt da R n u R* n (x)u* (x) = — I x . (79) 

J — oo J n= QJ-T J 1 — \R\ 

We can now combine this with Eq. (|50|) to write 

(V) = (hT2) 2(V)l ' (80) 

This result is the analog of Eq. ([71]) for the case of a Fabry-Perot cavity. In both cases, the 
momentum fluctuations grow as the square of the effective number of bounces. 

D. The Equal Arm Interferometer 

In an equal arm interferometer with a perfect (loss-free) 50-50 beam splitter, the input 
power is divided equally between the two arms. In the late 1970's, there was a controversy 
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over whether radiation pressure fluctuations will create noise in such an interferometer. The 
arguments reviewed at the beginning of this section leading to the standard quantum limit, 
Eq. (pop, assume that the radiation pressure fluctuations in the two arms are uncorrected. 
However, one would expect that a fluctuation which sends more power into one arm will 
cause a corresponding deficit in the other arm. This would lead to anticorrelated pressure 
fluctuations. Caves resolved this controversy in the context of the photon number 
approach. He showed that when vacuum modes which enter an unused port of the interfer- 
ometer are included, the fluctuations are uncorrelated. In this subsection, we will rederive 
this result using the stress tensor approach. 

1. Properties of a Beam Splitter 

Region 1 



u 




v 

FIG. 3. u and v represent the two incident wavepacket. r , t are the complex amplitude 
reflectivity and transmissivity for incident light u, and r' , t' are for the incident light v 

Let us consider a perfect beam splitter. (See Fig. It needs to satisfy the reciprocity 
relations, originally derived by Stokes in 1849, 

\r\ = \r'\ , \t\ = \t'\ , (81) 

M 2 +i*r=i (82) 

and 
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r'f + r*t' = 0, 



(83) 



where r , t are the complex amplitude reflectivity and transmissivity for light incident from 
one side and r' ,t' for light from the other side. The first equation, Eq. (|8"T|), arises form 
the assumption that the reflectivity is the same from both sides. The other two equations 
are the results due to the additional assumption of a no-loss beam splitter. Assume that 
we send in a wavepacket uq into the beam splitter in the x direction. Then it will reflect 
the amount tuq to region 1 and transmit tu$ to region 2. If this is a no-loss beam splitter, 
then the inverse operation will bring the reflected and transmitted wavepackets back to the 
original incident wavepacket, 

(ru o y + (tu )t* = u , (84) 

which leads to Eq. fl8"2"|). If there is another wavepacket U\ coming into the beam splitter 
from the y direction as well as the one Uq from x, then in region 1 we have the transmitted 
wavepacket i! U\ in additional to the reflected wavepacket ru , Similarly we get tu + r' u\ 
in the region 2. Again the inverse operation yields the relations 

(r m + t' «i)r* + (tu + r' u x )t* = u (85) 

and 



(r u + t' Mi)t'* + (t u + r' ui)r'* = u x (86) 

for a no-loss, perfect beam splitter. Both equations here will lead to the reciprocity relation 
Eq. (^). For a 50 — 50 beam splitter, r' = r, t' = t. Without losing generality, we may take 
the coefficients to be 

r = -^=e i4,r 

and 

t = -=e i<h > 

where <ft r , (fit are the phase change due to reflection and transmission, respectively. Plug 
them into the reciprocity relation Eq. (p3|), and the phase difference becomes 



TiTT 

A = <j) r -<p t = — ,n= 1,3,5... . (87) 

This phase difference is crucial in the discussion of the momentum correlation of the two 
end mirrors in interferometer. 



2. Correlation between two arms 

In a quantum mechanical treatment of the radiation field in the presence of a beam 
splitter, one often speaks of light entering one port of the interferometer and of vacuum 
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entering the other, unused port |l|,p|J10[| . From our point of view, this language is misleading. 
Vacuum modes are everywhere, and are entering both ports. Furthermore, there are an 
infinite number of such vacuum modes, and one would like to see more clearly which ones 
are actually relevant in a given situation. 

Assume that wavepacket uq (the occupied mode) with a particular frequency is incident 
on the beam splitter. It will reflect ruo to mirror #1 and transmit tuo to mirror #2 (See 
Fig. ^). Similarly, at mirror #1 there are modes r Uj reflected from the vacuum fields coming 
from the input port and tuu transmitted from the vacuum fields coming from the output 
port. At mirror #2, we have vacuum modes tuj and ruk in addition to the occupied mode, 
tuo- If we consider the momentum difference transferred to the mirrors p = Pi — P2, then 
the deviation becomes 



(Ap 2 ) = (Ap\) + (Apl) - 2(A Pl Ap 2 ) 



(88) 



Now in arm #1 the incident wavepacket is tuq and the complete set of wavepackets from 
vacuum are (J2j ru j + J2kt u k)- Follow the reasoning leading to Eq. (|5Q) in the one arm 
case; the momentum dispersion of mirror 7^1 becomes 

(A P 2 ) 
1 1 2 



R e {J dt J dar uo(^^r Uj + y ^tuk)* J dt' J da' T*UQ^^TUj + y^titfc)] 
da uqUq 



\z\ 
~2 



dt 
dt 



3 

2 1 



* _ J Re[(l + e jA )(l + e- 



)] 



da UqUq 



(89) 



Here we used the result A = ^ from Eq. (0). The momentum dispersion of mirror #2, 
(Ap|) will be the same as (Ap 2 ). The correlation between mirrors is 



(A Pl Ap 2 ) 



2 R- e lJ dt J daruo(^2rUj + ^2tUk)*J dt' J da't*ul(Y^t 



U ; 



k 



dt 



dau u* \ 2 * -Re[(l + e lA )(l + e lA )] 



\z\ 
0. 



J dt J dauQU* ^ * -Re[±2i] 



(90) 



We can see that there is no correlation between arms. The fluctuations are totally indepen- 
dent to each other. The dispersion of the momentum difference becomes 

(Ap 2 ) = 2(Apl) 
For b bounces, the dispersion is 



dt J dauQU*^j = 4(n)cu z 



(91) 



[Ap 2 ) b = 2b 2 (ApD =4b 2 (n)u 2 



(92) 



This confirms that Eq. 
interferometer. 



I]) give the correct velocity dispersion of each end mirror in the 
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V. DISCUSSION AND CONCLUSIONS 



In this paper, we have shown how quantum fluctuations of radiation pressure arise from 
fluctuations of the stress tensor operator. Our results are in agreement with those obtained 
previously using a photon number counting approach. In our approach, the radiation pres- 
sure fluctuations in a coherent state are due entirely to the cross term in the product of 
stress tensors. This term is both dependent upon the quantum state, and is singular in the 
limit of coincident points, however, we found that careful treatment of the integrals over 
space and time leads to a finite result. 

The cross term can be interpreted as representing the interference between vacuum fluc- 
tuations and the real photons present. Thus radiation pressure fluctuations in the stress 
tensor approach are driven by vacuum fluctuations. It is useful to compare the photon num- 
ber and stress tensor approaches at this point. Both approaches yield the same answers for 
all of the questions which were posed in this paper. (A possible exception is the Fabry- Perot 
cavity discussed in Sect. |1V C 2| using only the stress tensor approach.) However, the con- 
ceptual pictures presented by the two approaches are quite different. In the photon number 
approach, the pressure fluctuations on a single mirror are attributed to statistical variations 
in the numbers of photons striking the mirror. However, when one wants to treat the prob- 
lem of noise in an intererometer, especially the lack of correlation between the fluctuations 
in the two arms, it is necessary to invoke vacuum fluctuations Jl|||. In our view, the stress 
tensor approach provides a more unified description in which the role of vacuum fluctuations 
is clear from the outset. It is also likely to generalize more easily to complex situations. For 
example, all of the treatments of radiation pressure fluctuations, with which we are aware, 
assume that the end mirrors are perfectly reflecting. However, the stress tensor approach 
could be easily adapted to account for the finite reflectivity of this mirror. 

Radiation pressure fluctuations will play a role in laser interferometer detectors of gravity 
waves, especially in the future. At that point, it should become possible to measure these 
fluctuations experimentally. Confirmation of their existence can be viewed as experimental 
evidence for the reality of the cross term. 

Of special significance is the role of radiation pressure fluctuations in understanding the 
fundamental physics of stress tensor fluctuations. It seems natural that the same principles 
which apply to the stress tensor as a source of pressure on a mirror should also apply to the 
stress tensor as a source of gravity. We have seen that the cross term is essential to under- 
stand radiation pressure fluctuations. It then follows that the cross term must be included 
in the treatment of spacetime metric fluctuations driven by stress tensor fluctuations. 

Acknowledgement: We would like to thank K. Olum and C. M. Caves for valuable 
discussions. This work was supported in part by the National Science Foundation under 
Grant PHY-9800965. 

APPENDIX 

In this appendix, we calculate the integral J defined in Eq. fl34l) . Define u = t\ — 1% and 
v = ti + 1 2 . Next use the identities 
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and 



cos(c<jt 1 ) cos(ujt 2 ) = ~[cos(uju) + cos(uv)] (Al) 



dtidt 2 = - ( / du dv+ du dv ) . (A2) 

-T — w JO Ju 



Jo Jo 2 
After evaluation of the f-integrations, we may write 

r (r-u)(u 2 + a) 1 r u 2 + a r . . . n , . n . 

J= / aw — ^ — cosuju—— du-—. — sin urn + sinw n — 2r) . (A3) 

Jo (m 2 — o 2 )' 4 2cj jo (m 2 — o 2 ) J 

The second integral in the above expression approaches a constant as r — > oo, whereas the 
first integral contributes a linearly growing term: 

This integral contains third-order poles at u = ±b. It can be expressed as 

J = \ r{J+ + J-) , (A5) 

where in J + we assume Imfe > and in I + we take Im& > 0. Each of these integrals is in 
turn expressed as 

/± = ^(/±i + /± 2 ), (A6) 



where 



'±i 



00 , u 2 + a 



du 2 e- (A7) 



and 



'±2 



'OO 



oo 



du w^f e ~ wu - (A8) 



Each of these integrals is evaluated by closing the contour of integration in the appropriate 
half plane, and then evaluating the integral by a combination of integration by parts and 
Cauchy's theorem. For example, in the case of i+i, we close in the upper half plane and 
write 



2i-oo V rfM 2 u-bj (u + b) 3 2i-oo w- 



it — b du 2 



u 2 + a e<Wtt 



(u + o) ; 



= 717 



dw 2 



+ ° e ic« 



(u + by 



(A9) 



Similarly, we find 
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'+2 



— Til 



du 2 



u + a 
{u-bf 



(A10) 



I-i 



-Til 



du 2 



u 2 + a 

(u + by 



(All) 



u=b 



and 



1-2 = vri < — - 



{- 


u 2 + a iuJU 


) 


\du 2 


_(u-b) 3 





(A12) 



We may now combine all of these results to obtain the expression for J, Eq. (^). Note 
that this calculation involves integrations by parts very similar to those used in Ref. || and 
illustrated in Sect. |I| of the present paper. 
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